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The inter-layer direct Coulomb coupling is analyzed in a charge ordered superlattice bilayer model 
in which pairing is supported by inter-layer Coulomb energy gain (potential energy driven super- 
conductivity). The 2D pair-condensate can be characterized by a charge ordered state with a 
"checkerboard" like pattern seen by scanning tunneling microscopy. The 2D <S> 3D quantum phase 
transition of the hole-content at T c , supported by c-axis optical measurements, is also studied. The 
pair condensation might lead to the sharp decrease of the normal state c-axis anisotropy of the hole 
content and hence to the decrease of inter-layer dielectric screening. The drop of the c-axis dielectric 
screening can be the primary source of the condensation energy below T c . We find that a net gain in 
the electrostatic energy occurs along the c-axis, which is proportional to the measured condensation 
energy (Uo) and with T c : E^ D » 2(£ a f,/ao + l) 2 ?7o ~ ksT c and is due to inter-layer charge com- 
plementarity (charge asymmetry of the boson condensate) where ^ a b is the coherence length of the 
condensate and ao ~ 3.9A is the in-plane lattice constant. The model naturally leads to the effective 
mass of m* ~ 4m e found by experiment. The static c-axis dielectric constant e c is calculated for 
various cuprates and compared with the available experimental data. We find correlation between 
T c and the inter-layer spacing d, e c and with the coherence area of the condensate. 
PACS numbers: 74.20.-z, 74.25.-q, 74.72. -h 
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I. INTRODUCTION 

It is more or less generally accepted now that the con- 
ventional electron-phonon pairing mechanism cannot ex- 
plain cuprate superconductivity, because as high a tran- 
sition temperature as 16AK (the record T c up to now 
[1]) cannot be explained by the energy scale of lattice 
vibrations without leading to lattice instability [2]. It 
is already well established that much of the physics re- 
lated to high-temperature superconductivity (HTSC) is 
in 2D nature, one of the basic questions to be answered, 
however, in the future is whether HTSC is a strictly 2D 
phenomenon or should also be described by a 3D theory. 

A well known experimental fact is that the zero re- 
sistance occurs along the ab plane and the c-axis at the 
same critical temperature [3] suggesting that there must 
be inter-layer (IL) coupling involved in the mechanism 
which drives the system into HTSC. Another experiment, 
such as the intercalation on Bi2212 [4], however leads 
to the opposite conclusion. Intercalation of I or organic 
molecules, which expands the unit cell significantly along 
the c-axis, does not affect T c . This finding is against 
IL coupling and supports low-dimensional theories. The 
large anisotropy of the resistivity (and of other transport 
properties) is again not in favour of 3D theories of HTSC 
[3]. 

Experiments on Y 'Ba-iCu^O^ - y (YBCO) ultrathin ar- 
tificial HTSC compounds, sandwiched between thick 
nonsuperconducting PrBa 2 Cu30r- y (PBCO) layers 
[3,5] and measurements on (BaCu0 2 + x ) 2 / '(CaCu0 2 ) n 
heterostructures [6] indicate the continous decrease of 



T c with the decreasing number of superconducting 
(SC) layers. In particular the one-unit-cell thick sam- 
ple of YBCO/PBCO heterostructure exhibits T c w 
20 K [5,7]. Other heterostructures, such as the 
(Bao.9Ndo.iCu0 2+x ) 5 /(CaCu0 2 )2/(Ba . 9 Nd .- i Cu0 2+x ) 5 , 
containing a single bilayer SC block isolated from each 
other by insulating blocks, were shown to have T c w 55 
K [8]. These findings again indicate the importance of 
IL coupling in high temperature superconductors. There 
are a couple of other findings against and supporting the 
3D nature of HTSC [9-11]. Most notably the systematic 
dependence of the transition temperature T c on the o 
axis structure and, in particular, on the number of Cu0 2 
planes in multilayer blocks are also strongly in favour of 
the 3D character of HTSC. It is therefore, a fundamental 
question whether at least a weak IL coupling is needed 
for driving the system to a superconducting (SC) state 
or a single Cu0 2 layer is sufficient for HTSC [12]. 

There has been considerable effort spent on under- 
standing HTSC within the context of IL coupling mecha- 
nism in the last decades in which c-axis energy is available 
as a pairing mechanism [13-16]. In other approaches the 
importance of IL hopping is emphasized vs. the direct 
IL Coulomb interaction of charged sheets [9,12]. There 
is another theory providing explanation for HTSC using 
the general framework of BCS combined with IL coupling 
[17]. The so-called IL tunnelling (ILT) theory [9,12], how- 
ever, is no longer considered a viable mechanism for SC 
in cuprates since ILT could provide no more then 1% of 
the condensation energy in certain cuprates [2,10,18]. 

In this paper we propose a simple phenomcnologi- 
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cal model for explaining the 3D character of HTSC in 
cuprates supported by calculations. We would like to 
study the magnitude of direct Coulomb interaction be- 
tween charge ordered square superlattice layers as a pos- 
sible source of pairing interaction. Our intention is to 
understand HTSC within the context of an IL Coulomb- 
mediated mechanism. Of particular relevance to our in- 
vestigation are the doping, multilayer, and pressure de- 
pendence of T c in terms of a charge ordered superlattice 
nature of pair condensation. 

The size of a charge ordered characteristic superlattice 
can directly be related to the in-plane coherence length 
£ a b of cuprates, which is proportional to the linear size of 
the pair condensate (real space pairing) in the ab-plane 
[19]. The relevant length scale for superconductors is the 
characteristic size of the Cooper-pair and can be esti- 
mated by un uncertainty- principle argument [19]. Only 
those charge carriers play a decisive role in HTSC which 
has energy within <~ ksT c of the Fermi energy and sets in 
at T c with the momentum range Ap w ksTc/vp, where 
vp is the Fermi velocity, leading to the definition of the 
Pippard's characteristic length [19] 



cuprates can be summarized in the parabolic dependence 
of the maximum critical temperature T™ , 



T c /T c m -l-82.6(p-p ) 2 , 



(2) 



TlVF 

'k B T c 



(1) 



where a is a numerical constant of order unity, to be 
determined. £ a (, is a relevant number for describing the 
2D confinement of the pair-condensate wave function be- 
low T c . £ a b can also directly be obtained according to 
the anisotropic Ginsburg-Landau theory via the measure- 
ment of the upper critical field £f C 2,a& [19]. It has been 
shown recently that £ a f> obtained from Eq. (1) or from 
H C 2,ab measurements are very close to each other [20]. 

The IL charging energy we wish to calculate depends 
then on the IL spacing (d), the IL dielectric constant 
e c , the hole content p and the size of the superlattice. 
A large body of experimental data are collected which 
support our model. A consistent picture is emerged on 
the basis of the careful analysis of this data set. Finally 
we calculate the static c-axis dielectric constant e c for 
various cuprates which are compared with the available 
experimental observations. 



II. OPTIMAL DOPING AND THE 2D-3D PHASE 
TRANSITION OF THE HOLE-CONTENT 

It is commonly accepted that charge carriers are 
mainly confined to the 2D CuOi layers and their con- 
centration is strongly influenced by the doping agent via 
hole doping. Holes (no charge and spin at a lattice site) 
in the 2D CuOi layers are the key superconducting el- 
ements in high temperature superconductivity (HTSC). 
A characteristic feature of many high temperature super- 
conductors (HTSCs) is the optimal hole content value of 
p w 0.16e per CuOi layer at optimal doping measured 
in the normal state (NS) [21]. This general feature of 



where T™ corresponds to the optimal hole concentration 
Po [22]. Since T c appears to be maximized at p w 0.16, 
we pay special attention to the IL Coulomb interaction 
in cuprates at this particular hole concentration. 

Our starting hypothesis is that the hole content goes 
through a reversible 2D <^ 3D quantum phase transition 
at T c in layered cupper oxides. Optical studies on the c- 
axis charge dynamics reveals this phenomenon: the c-axis 
reflectance is nearly insulating in the NS but below T c is 
dominated by the Josephson-like plasma edge [23] . Below 
T c a sharp reflectivity edge is found at very low frequen- 
cies (lower then the superconducting gap) for a variety 
of cuprates arising from the carriers condensed in the SC 
state to the afr-plane and due to the onset of a coherent 
charge transport along the c-axis [24-28]. The appear- 
ance of plasma in the SC state and the absence of it in 
the NS seem to be a common feature of HTSC cuprates 
[27]. Band structure calculations predict an appreciable 
c-axis dispersion of bands close to the Fermi surface and 
thus an anisotropic three-dimensional metallic state [29] . 
Other first principles calculations indicate that above T c 
the hole charge p is charge transferred to the doping site 
in the charge reservoir (2D — > 3D transition) [30,31]. 

Furthermore, we expect that below T c the hole-content 
is condensed to the sheets forming anti-hole regions (hole- 
content charge at a lattice site, 3D — ► 2D transition) . An 
anti-hole corresponds to an excess charge condensed to a 
hole lattice site in the sheets below T c . Therefore, in the 
normal state hole doping is the dominant charge transfer 
mechanism while in the superconducting (SC) state the 
opposite is true: the IL hole-charge is transferred back 
to the CuOi planes (charging of the sheets). The latter 
mechanism can be seen by the measurement of transport 
properties as a function of temperature [3,32-34] if we 
assume that the increase in the density of in-plane free 
carriers below T c is due to pair condensation. The sharp 
temperature dependence of the c-axis dielectric constant 
e c and optical conductivity [35] seen in many cuprates 
and in other perovskite materials also raises the possibil- 
ity of a 2D 3D charge density condensation mecha- 
nism at T c [3,38]. Therefore, the pair condensation can be 
described by an anisotropic 3D condensation mechanism, 
and by a doping induced 2D-3D dimensional crossover 
[11]. ■ 

In any naive model of electron pairing in cuprates the 
Coulomb repulsion is troublesome. When the pairs of 
charged carriers are confined to the sheets, naturally 
a net self-repulsion of the pair condensate occurs. Al- 
though short-range Coulomb screening in the dielectric 
crystal can reduce the magnitude of the repulsion but it 
is insufficient to cancel Coulomb repulsion [39]. Spatial 
separation should reduce the interaction strength, but 
even at 1AA e 2 /r w leV if unscreened. One can as- 
sume capacitative effect between the Cu0 2 planes: The 
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charged boson condensate in one plane is stabilized by 
a deficiency of that charge on another plane [39]. The 
various forms of the capacitor model are associated with 
the 3D character of HTSC, considering the inter-layer 
charge reservoir as a dielectric medium [39]. The inter- 
layer charging energy might be insufficient to stabilize 
the self-repulsion of the holes in one plane and the self- 
repulsion of the charge condensate on another plane in 
the capacitor model. Also, the superconducting proper- 
ties of the hole-rich and nearly hole-free sheets would be 
different which is not verified by experimental techniques. 

Our intention is to combine the 2D and 3D nature of 
HTSC. Therefore, we propose a phcnomcnological model 
in which the charge distribution of the planes is polar- 
ized in such a way that holes and anti-holes (hole-electron 
pair) are phase separated within each of the sheets lead- 
ing to a charge ordered state (COS). Recently, signatures 
of charge ordering have been found in various cuprates 
and manganites in the presence and the absence of mag- 
netic fields [41-44]. Furthermore, the "checkerboard" 
charge pattern seen in Bi2212 [41] with spin periodicity 
of ~ 8a strengthens our expectation that the SC state 
can be characterized by a COS. We focus on, then the 
magnitude of IL coupling (direct IL Coulomb interaction) 
between 2D static charge ordered superlattices. 

Charging of the sheets: The CuOi layers carry neg- 
ative charge obtained from the charge reservoir even in 
the insulating stochiometric materials, e.g. in the infi- 
nite layer compound CaCuOi each layer is charged by 
2e charge donated by the Ca atoms. In (Cu02) 2 ~ then 
the electron configuration of Cu and O is 3(i 10 4s 1 and 
2p 2 z . Unit (CuC>2) 2 ~ is typical of any undoped (stochio- 
metric) cuprates and ab initio calculations provide ap- 
proximately the charge state Cu 15+ C>2 ' 5 ~ [3,29] which 
is due to charge redistribution between Cu and O. The 
{CuO-i) 2 ^ plane itself is antiferromagnetic and insulat- 
ing, e.g. there are no holes in the "overcharged" layers. 
Upon e.g. oxygen doping, however, the doping charge p 
is transferred along the c-axis to the doping site in the 
charge reservoir, since the doping atom exhibits a rare 
gas electron structure (0 2 ~, if oxygen doping) [30,31]. 

Briefly, the c-axis anisotropy of the hole content is 
strongly temperature dependent in cuprates. Our ba- 
sic assumption is that in the SC state every second hole 
is filled up with the charge of 2p due to the condensation 
of the hole-content to the sheets leading to a charge or- 
dered state in the SC state. Therefore the annihilation 
of the half of the holes results in the appearance of a 
hole-anti-hole pair. The hole-anti-hole pairs are the ba- 
sic quasiparticles of the SC COS. The phase separation 
of the holes and anti-holes (leading to a charge ordered 
state) is stabilized by the intra- and inter-layer hole-anti- 
hole interactions. Basically the phase separation of holes 
and anti-holes in the planes is the manifestation of strong 
Coulomb correlation in the SC state. The supercon- 
ducting charge-ordered states in cuprates are character- 
ized theoretically based on the assumption that cuprates 
are in close proximity to quantum critical points where 
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FIG. 1. Upper panels: The charge odered state of the type 
of a "checkerboard" of the hole-anti- hole condensate. Opened 
and filled circles denote the holes with a charge of qi — +0.16e 
and anti-holes (qj = — 0.16e), respectively in the 4ao x 4ao 
(5 x 5) square lattice layer model. Note that the left panel 
accomodates 13 anti-holes which corresponds to 2e + 0.16e/2 
charge. The right panel contains 12 anti-holes corresponding 
to 2e — 0.16e/2 charge. The two lower panels correspond to the 
antiferromagnetic insulating state (left) and to a non charge 
ordered holed-doped system (normal state, right). 




FIG. 2. The charge ordered state of the hole-anti-hole con- 
densate in the 4ao x 4ao charge ordered bilayer superlattice 
model. Note the charge asymmetry between the adjacent lay- 
ers. The bilayer can accomodate a pair of boson condensate 
(4e). Noteworthy that holes (empty circles) and anti- holes 
(filled circles) can be characterized by stripes along the diag- 
onal lines. The inter-layer charge complementarity of these 
charge stripes is crucial for the inter-layer Coulomb energy 
gain. 
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spin/or static charge order occurs in the SC state [45]. 
In this paper we will focus on the Coulomb interaction 
between the charge ordered sheets screened by the IL di- 
electric media (charge reservoir). It will be immediately 
apparent from our analysis that net energy gain occurs 
due to the IL Coulomb interaction of the charge ordered 
layers when the hole-electron charge pattern is asymmet- 
rically distributed in the adjacent layers (Fig. 2). 

III. THE SUPERLATTICE MODEL: THE 
COHERENCE AREA 

We propose to examine the following superlattice 
model of pair condensation: A pair of charge carriers 
(2e) van be distributed over 2/0.16 = 12.5 CuOi unit 
cells in a square lattice layer if the 2e pair is composed 
of the hole content po w O.IQ/C11O2 at optimal doping. 
However, allowing the phase separation of hole-anti-hole 
pairs, every second unit cell is occupied by — 0.16e (anti- 
hole), and the rest is empty (holes, +0.16e), therefore we 
have 25 unit cells for a condensed pair of charge carriers 
(Fig 1., that is the unit lattice of the pair condensate 
in the 5x5 supercell model). Therefore, 2po — — 0.32e 
hole charge condenses to every second hole forming anti- 
hole sites. The 25 lattice sites provide the hole content 
of 25 x 0.16e = 4e and therefore 2e excess charge in the 
anti-hole sites. In other words a Cooper wave-function 
can be distributed on a 5 x 5 superlattice with a node of 
— 0.16e in every anti-hole site. The remarkable feature is 
that the size of the 5x5 condensate (four lattice spacings, 
4a w 15. 5A) is comparable with the measured small co- 
herence length of single-layer cuprates (£ a b ~ lOA to 
20A) [3,19,44,46]. £ a b can directly be related to the char- 
acteristic size of the wave-pocket of the local Cooper pair 
(coherence area) [3,19]. The charge ordered superlattice 
model can in principle be applied not only for po ~ 0.16e 
but also for the entire doping regime. The only difference 
what happens is that the charge nodes in the charge or- 
dered state changes upon doping and consequently the 
inter-layer coupling scales with the anti-hole charges in 
the underdoped regime. In the overdoped regime, how- 
ever, the situation is more difficult [51,57]. 

Our expectation is that the COS of the 5x5 model 
given in Fig. 1 can be an effective model state for de- 
scribing the SC state. An important feature of this model 
is the charge separation dq in the charge ordered state, 
where +0.16e and — 0.16e partial boson charges are local- 
ized alternatively (dq = 0.32e). The hopping of charge 
carriers from the anti-hole sites to the holes can reduce 
the magnitude of the charge separation dq leading to the 
extreme case when dq w which is nothing else then the 
(Cu02) 2 ~ antiferromagnetic insulating state. Therefore 
characteristic quantity of the SC state dq is directly re- 
lated to the hole-content seen in cuprates in the NS as a 
function of doping. 

The hopping of anti-holes in such a charge ordered lat- 
tice layer might lead to SC without considering any lat- 



tice vibration effects when dq is in the optimal regime. 
The freezing of the COS leads to insulating Wigner crys- 
tal (the NS) [40]. Therefore the melting of the crystal- 
lized COS is required for SC [45]. HTSC can be char- 
acterized in this way as a competition between a liquid- 
crystal-like striped COS (SC) and a Wigner-crystal-likc 
phase (NS) in a semiclassical theory of hole dynamics 
[40]. 

Below T c the charged-ordered state becomes stable 
compared with the competing phase of the NS supported 
by IL charging energy. Holes and anti-holes are placed in 
such a way in adjacent layers to maximize IL charging en- 
ergy. Therefore, holes in one of the layers are always cov- 
ered by anti- holes in the other layer and vice versa (FIG. 
2, inter-layer electrostatic complementarity, bilayer 5x5 
(4a x 4a ) model). An important feature is then that 
the boson condensate can be described by an inter-layer 
charge asymmetry. The IL coupling of the boson-boson 
pairs in the bilayer 5x5 model naturally suggests the 
effective mass of charge carriers m* ss 4m e , as it was 
found by measurements [47,50]. In the next section we 
generalize the 5x5 model to represent a real space pe- 
riodicity of N x N coherence area. The charge ordered 
state presented in this article can also be studied as a 
charge density wave (CDW) condensed on a superlattice. 
Using the CDW terminology one has approximately two 
CDWs within a characteristic bilayer with an amplitude 
of <~ 0.16e. This system might not be insulating if we 
assume high rate of intersite hopping and low rate of 
IL hopping of the charge carriers. Then the charge or- 
dered state is stabilized due to the attractive IL Coulomb 
interactions of asymmetrically condensed neighbouring 
CDWs. 

We characterize the purely hole-doped NS as follows: 
In each sheet each CuO-i units posses p = +0.16e charge. 
The hole charge p is charge transfered to the charge reser- 
voir (doping site) above T c . In this model one can naively 
expect that the self-repulsion of a hole rich sheet can be 
quite large. In order to reduce the hole-hole self-repulsion 
of the sheets, a hopping of ~ p charge between the ad- 
jacent hole sites would result in the reduction of the re- 
pulsion leading to a Wigner-crystal-like COS depicted on 
FIG 3. As is well known, the hopping of certain amount 
of charge (~ p) between the holes leads to metallic con- 
ductance which is typical of hole-doped cuprates. This is 
the normal state (NS) of the cuprates. Due to the strong 
c-axis anisotropy of the hole content in the NS IL dielec- 
tric screening is large in this state because of the increase 
in the dynamic component of the IL dielectric permittiv- 
ity (e c ) and therefore the IL interaction energy vanishes, 
E^ s ss 0. That is the enhanced dielectric screening in 
the NS due to the localization of the hole content in the 
charge reservoir. The lack of the IL plasma edge in the 
c-axis optical spectra of various cuprates in the NS is a 
strong evidence of this phenomenon [60] . 
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FIG. 3. A possible charge ordered state of the normal state. 
Opened and filled circles denote the holes with a charge of 
qi « +0.32e and anti-holes (qj w — 0.16e), respectively in the 
5x5 square lattice layer model. 

Our proposal is to understand HTSC in such a way 
that the pair-condensate is stabilized at higher tempera- 
ture by the IL-charging energy provided by the antisym- 
metry of the pair-condensate charge density between the 
adjacent layers. The c-axis dielectric constant of the SC 
state is reduced to the average value of the IL dielectric 
medium due to the pair condensation of the hole charge 
to the sheets. No IL screening of the hole charge occurs 
in this state which may well lead to energy gain . 



IV. THE TOTAL ENERGY OF THE PAIR 
CONDENSATE 

Our intention is to develop a simple working hypoth- 
esis in which the IL charge reservoir provides an av- 
erage dielectric background and the hole content fur- 
ther enhances IL dielectric screening when 2D — > 3D 
transition occurs (hole-doping, the NS). The dielectric 
plasma provided by the hole charge results in the dy- 
namic screening effect of Coulomb interaction which is 
typical of hole doped cuprates in the NS. In the oppo- 
site case (3D — > 2D, pair condensation) the IL dielectric 
screening is nearly reduced to the average background 
value of the IL ion core spacer. This is what leads to 
IL energy gain. The underlying source of the condensa- 
tion energy is then the energy gain due to the lack of 
dynamic screening in the SC state. The possibility of di- 
rect IL hopping of Cooper pairs is not considered within 
this model as it was ruled out as important aspect of 
HTSC [10]. 

We start from a very general description of our model 
system using e.g. a Hamiltonian similar to that is given 
elsewhere [14] or which can also partly be seen in gen- 
eral text books [58]. We would like to describe then the 
2D 3D condensation of the hole-content using the 
Hamiltonian 



(3) 



where Hf D is the BCS-type Hamiltonian of the intra- 
layer condensate. 



fl? D = £ c * c L<w (4) 

k,a 

+ V J2 C U,i C -kl,i c kl,i c -k1,i- 
k,k' 

In Eq. (4) cj^ i is the creation operator for electrons 
in the ith layer, with linear momentum k within the 
layer and spin a, using the effective- mass approximation 
€k = h 2 k 2 /2m*. V = —\V\ is the attractive inter-layer in- 
teraction (3D coupling) and is assumed to originate from 
some of the proposed mechanisms [9,13-15]. The attrac- 
tive IL term HfP given by 

H if = ~ t 4a,i C k^+Yj2 C L,i C -k/3,j C -kl3,jCka,i (5) 

k(j,ot : fi k.k' 

+ W £ C la,i C -k0,i C -k0,jCk a ,j- 
k,k',a,/3 

The first term in Eq. (5) is the direct IL hopping, while 
the second and third terms describe IL coupling in a 
general form, t is very small in oxide superconductors 
[10,14]. Y denotes the IL coupling assisted by direct 
Coulomb interaction between charged layers (this is of 
particular interest in our model). W is the coupling 
constant and can arise through the Coulomb interac- 
tion causing inter-band transitions at the Fermi surface 
to some of the occupied or empty bands away from the 
Fermi level, with finite dispersion, however, along the c- 
axis. In this paper we consider only the second term in 
Eq. (5) as the source of the attractive interaction for 
pairing and neglect the rest of the inter-layer Hamilto- 
nian HfP (t«0,WwO). 

We calculate then the energy of the nearly 2D electron 
pair condensate in the CuOi plane for the bilayer super- 
cell problem. It is assumed that the pair condensate be- 
haves as a nearly free electron gas with afe-kinetic energy 
Tfo n and potential energy, which is mainly its in-plane 
self-electrostatic energy (E^ b ) and the out-of-plane inter- 
layer interaction energy E* L . For simplicity, the rest of 
the electron system is neglected completely. The exter- 
nal potential of the condensate is also excluded in this 
model, that is the lattice-condensate interaction, which 
is assumed to be negligible in the charged (CuC>2) 2 ~ sys- 
tem (at least its contribution is negligible to the conden- 
sation energy). In other words the ionic background of 
the planes is screened by the core and valence electrons 
of the (CuC>2) 2 ~ plane. The kinetic energy of the charge 
condensate is due to the hopping of the charge carriers 
between the CuOi sites within the sheets. We do not take 
into account the complications due to ionic heterogene- 
ity, nonpointlike polarization, etc. The effective Coulomb 
interaction between spinless point-charges may be ap- 
proximated by the expression V e //(r) = e 2 /(47re e c r), 
where e c takes into account phenomenologically the di- 
electric screening effect of the IL dielectric medium and 
confined hole charge. This kind of a rough approxima- 
tion has widely used by several groups in the last decade 
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[13,39,47,70]. The important feature is that the c-axis 
dielectric screening (e c ) is nearly reduced to the static 
value of the average background dielectric constant in 
the SC state. In accordance with this a sharp temper- 
ature dependence of e c (u) is found in BSCCO by c-axis 
optical measurements at T c [35]. The static out-of-plane 
dielectric function can be obtained from the sum ride 



7T Jn OJ' 



(6) 



where £2(0;) is the dynamic component of e c (u>) [68]. 

Although, the r-dependence of the screened in-plane 
electrostatic interaction between the condensed charges 
Qi and qj is not perfectly 1/r, we approximate it with 
the expression V£fj ~ e 2 /47reoe a ;,r as well, where screen- 
ing is taken into account implicitly via e a b- The lowest 
eigenvalue of H given in Eq. (3) is E = (\I>|.ff|\I>), where 

^{fi, T2) — ^2^g{k)e tk ' ri e lk ' T2 being the wavefunction of 
the interacting pair and g(k) is a pair-correlation func- 
tion [19]. Then the energy of the lowest eigenstate of the 
condensate in the SC state is approximated by using a 
pure dielectric form for the potential energy 



E,, 



m+Ef + E? 



■A*(n,r 2 ) 



+ 



2m* 

fl 2^ (l) 
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1 2 00 2N 2 (n) (to) 



n=l m—2 ij 



(n,m) 



where % is the Planck constant, m* ~ 4m e is the effec- 
tive mass for the holes induced in the half-filled bands 
[47], n, m are layer indices, and rj"' m ^ are the intra- 
layer and inter-layer point charge distances, respectively. 

The most important components are the interactions 
( 1 2) 

with r\j' (bilayer components), however one has to 

sum up for the IL interactions with terms r|™'"^ , where 
m = [2,oo]. Note that only the interactions of various 
layers with the basel bilayer (FIG. 2) are considered along 
the c-axis in both directions (up and down). The in-plane 
electrostatic screening e a b is completely separated from 
the out-of-plane dielectric screening (e c ). qi,qj are the 
partial point charges/atoms in the N x N superlattice 
model at optimal doping. 



= ± 



4e 



(8) 



The factor 3 is included in the denominator because each 
CuOi site consists of 3 atoms (qij w ±0.053e at optimal 
charge separation, N = 5). For the sake of simplicity it is 
assumed that the charges are equally distributed among 
Cu and O atoms within a CuC>2 site. The number of the 
lattice sites in the characteristic superlattice is N 2 and 



U « {N - l)oo 



(9) 



where ao w 3. 88 A is the ab lattice constant. Note that 
the coherence length £ a f> can be given in a unit and the 
real space period N can be compared with £ a b via Eq. (9) 
(see also Fig. (1)). Nh — 2e is the charge of the electron 
pair. The anti-hole charges qf hole must satisfy the charge 
sume rule within a characteristic bilayer over a coherence 



N 

E 



q ahole = 



(10) 



where q? hole represents the anti-hole point charges. Nat- 



urally the charge neutrality Y^if 



-.hole 



iholc 



) 





is also required. In this paper we study the lattice size 
N x N = 5, which we found nearly optimal for a vari- 
ety of cuprates. However, it can be interesting to study 
the variation of the " characteristic" lattice size in differ- 
ent cuprates as a function of various parameters (dop- 
ing, pressure etc.). The kinetic energy of the boson 
condensate arises from the hopping of anti-holes (hole- 
anti-hole exchange) between the adjacent sites within the 
sheets, against the electrostatic background of the rest of 
the hole-anti-hole system. Single hole-anti-hole exchange 
is forbidden since extraordinary repulsion occurs which 
leads to the redistribution of the entire hole-anti-hole 
charge pattern. The collective intersite anti-hole hopping 
results in the kinetic energy of the condensate. 



V. THE CONDENSATION ENERGY IN THE 
CHARGE ORDERED BILAYER SUPERLATTICE 
MODEL 

Important consequence of the model outlined in this 
paper is that the condensation energy (Uo) of the SC 
state can be calculated. We will show here that within 
our model the primary source of Uq is the net energy gain 
in IL Coulomb energy occurs due to the asymmetrical 
distribution of the condensed hole-charge in the adjacent 
layers (see Fig. 2) below T c . 

Per definition Uo is the free energy difference of the nor- 
mal state and the superconducting state [10]. Therefore, 
taking the energy difference U « AE tot = \E^J - E%£\, 
using Eq. (7) we get 



U w AE ab + AE l 



? IL 



(11) 



Tha a6-plane contribution to the SC condensation en- 
ergy AE ab is expected to be negligible at optimal dop- 
ing. That is because our Coulomb energy calculations 
indicate that the a6-plane Coulomb energy in the NS and 
in the SC state is nearly equal and therefore AE® b ~ 
and also the kinetic energy of the pair-condensate E% b n 
is much smaller by several order of magnitude then IL 
coupling. For instance we get for the prototypical 5x5 
lattice ~ —1.2 x 10 -20 J afe-Coulomb energy both for the 
NS and for the SC COS. The COS of the NS is the one 



G 



given in Fig (3). The calculated ab- kinetic energy of the 
5x5 codensate is ~ 1.1 x 1CT 29 J if we use the very 
simple formula (" electron in a box" ) 



- kin 



r 2e 



(12) 



where £ c is the c-axis coherence length. The inter-layer 
Coulomb interaction (attraction) E^ L w — 1.3 x 10~ 21 
J for the 5x5 bilayer. It must be emphasized, how- 
ever, that away from optimal doping the Coulomb en- 
ergy of the condensate might affect the magnitude of the 
condensation energy The abrupt jump of the measured 
condensation energy Uq seen in the slightly overdoped 
regime [56] can be attributed to the increased AE ab con- 
tribution to {Jo- 
in this article we restrict our study, however, to opti- 
mal doping and the study of the doping dependence of 
Uo is planned in the near future. The main contribution 
to Uo is then the IL energy gain AE[ L in the SC state 
can be given as follows, 



U « \E, 



IL,NS 



-E. 



IL,SC\ 



(13) 



We can further simplify Eq. (13) if the NS contribution 
to Eq. (13) is E£ L ' NS w 0, which holds if IL coupling is 
screened effectively in the NS (large density of the hole 
content in the IL space, large e c (NS)). Again, if we as- 
sume e c w 33, we get E£ L ' NS w +2 x 1(T 18 J (12.3 eV) 
for a charge ordered 5x5 lattice which would be an ex- 
traordinarily large value for the IL Coulomb repulsion. 
Assuming, however, e c w 10000, we get the more real- 
istic Coulomb repulsion of Ef. L > NS w +10~ 22 J. Indeed 
there are measurements for cuprates which indicate large 
e c in the NS in the range of 10 3 to 10 5 [63]. Extraordinar- 
ily large dielectric constant has also been found recently 
in perovskite materials [36]. Anyhow the "relaxation" of 
the huge IL Coulomb repulsion in hole-doped cuprates 
in the NS can not easily be understood without the con- 
sideration of a large e c . Nevertheless experimental mea- 
surements and theoretical speculations suggest that e c in 
cuprates and in other materials with perovskite structure 
is strongly temperature and doping dependent [35] . The 
sharp decrease of e c below T c should also be explained 
within our bilayer model by the pair-condensation of the 
hole-content to the sheets. This is what can be seen in 
the c-axis optical spectra of cuprates [23,60]. The lack 
of the c-axis plasma edge in the NS is an obvious exper- 
imental evidence of the strong temperature dependence 
of e r . Then we have 



2(n+l)A 2 [/ ~£, 



IL,SC 



(14) 



where Ef. L ' SC is the Coulomb energy gain in the SC state. 
Uo is the experimental condensation energy given per 
unit cell. Eq. (14) is generalized for multilayer cuprates 
introducing n, the number of Cu0 2 layers within a mul- 
tilayer block. 



The hole-conductivity in the SC state is strictly 2D 
phenomenon, no direct IL hopping of quasiparticles is 
considered within this model. T c is mainly determined 
by the inter-plane distance and by the static c-axis di- 
electric constant (the c-axis component of the dielectric 
tensor). An interesting feature of the NxN bilayer model 
is then that it is capable of retaining the 2D character of 
superconductivity while T c is enhanced by 3D Coulomb 
interactions. 



VI. THE CONDENSATION ENERGY AND T c 

There are number of evidences are available which 
suggest that HTSC occurs beyond the BCS limit. In 
those materials which contain nearly isolated single lay- 
ers, such as Bi 2 Ba2CuOe + s (T c = 20 K [69]) or su- 
perlattice structures (periodic artificially layered materi- 
als) such as O-doped {BaCuO^I (CaCuO<2) n thin films, 
when n = 1 [6] and in Y ' BaiCu^Oi-s / 1 PrBaiCu-sOi-?, 
(YBCO/PBCO) [3,5,7] the critical temperature is lim- 
ited to T c < 30K (below the BCS limit), therefore these 
materials are not considered as high-T c superconductors 
in this article. In these compounds the IL distance is so 
large that the Cu0 2 planes are decoupled and the super- 
conducting properties can be understood within the BCS 
formalism. However, the multilayer Bi-compounds (with 
the same charge reservoir as in the single layer Bi2201), 
the [BaCuO^ij {CaCuOi) n thin films, when n > 2 [6] 
or YBCO with thin PrBa 2 Cu 3 7 s layer [3,5,7] exhibit 
HTSC. Therefore it is worth to explain the enhancement 
of T c beyond the BCS limit assuming other mechanism 
than the elcctron-phonon coupling. Inter-layer Coulomb 
coupling can be a natural source of the condensation 
energy and of HTSC. Checkerboard-like charge pattern 
(COS) seen experimentally [41-43] directly leads to IL 
energy gain and to potential energy driven condensation 
energy if the hole-anti hole charge pattern is asymmc- 
trycally condensed to the adjacent layers (Fig 2). As- 
suming that thermical equilibrium occurs at T c for the 
competing charge ordered phases of the NS and the SC 
state the following equation for the condensation energy 
can be formulated using Eq. (14), 



2(n+l)A 2 [/ ~2(n+l) 



u + 1 

a 



U Q ^k B T c . (15) 



At a first look this formula seems to be unusual because of 
the dependence of the condensation energy on T c . The 
available measurements of the condensation energy on 
various cuprates show no correlation of Uo with the crit- 
ical temperature. 

One of the important goals of this paper, however, to 
show that correlation can indeed be found with T c if T c is 
plotted against 2N 2 Uo (the bilayer condensation energy). 
In other words the condensation energy of the coherence 
bilayer-hole system shows correlation with T c . In order 
to test the validity of Eq. (15) we estimate the coherence 
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TABLE I. The calculated coherence length of the pair con- 
densate using the experimental condensation energies of var- 
ious cuprates and Eq. (16) at optimal doping. 





T c (K) k B 


T c (meV) U (ueV/u.c 


) Cab C {ao) 




LSCO 


39 


2.5 


21 a 


~ 7 


5-8 h 


T12201 


85 


7 


100 ± 20 6 


~ 5 




Hgl201 


95 


7.8 


80 - 107 c 


~ 5 


5 C 


YBCO 


92 


7.5 


110 d 


~ 3 


3-4 9 


Bi2212 


89 


7.3 


107 e 


4 


4 - 5 f 



Uo is the measured condensation energy of various cuprates 
in /i eV per unit cell at optimal doping. a Uo ~ 2 J/mol 
from [52,53], b [18], c Uo « 12-16 mj/g from [54,67] 
and Z ab from [55], d U « 11 J/mol from [48,56,57], e 
Uo ~ 10 J/mol from [51], ' from recent measurements of 
Wang et al. [20], 9 from [19,49], ,l from [19], Cab" is cal- 
culated according to Eq. (16) and is also given in Table I 
and is the measured in-plane coherence length given 

in ao ~ 3.9A. The notations are as follows for the com- 
pounds: LSCO (Lai.ssSro.isGuOi), T12201 (ThBa 2 Cu0 6 ), 
Hgl201 (HgBa 2 CuO i+ s), YBCO {YBa 2 Cu0 7 ) and B12212 
is Bi2Sr2CaCu20s+s- 



length of the pair condensate using Eqs. (15) and (9) 
using only experimental data, 



£« b 



no 



2(n+l)Uo 



(16) 



The results are given in Table I as and compared 
with the available measured '. The agreement is ex- 
cellent which strongly suggests that Eq. (16) should also 
work for other cuprates. The validity of Eq. (15) is clearly 
clarified in FIG 4 using mostly experimental informa- 
tion for t; a b, T c and U$ (values are given in Table I.). 
For T12201 no measured t; a b is found in the literature, 
therefore the estimated t; a b is used (given in Table I.). 
In Section VII we will show that using this estimated 
value of N = 6 we predict e c in nice agreement with 
optical measurements for T12201. Remarkable feature of 
FIG 4 is that the slope of the linear fit to the measured 
data points precisely gives us ks which nicely confirms 
Eq. (15). The underlying physics of HTSC seems to be 
reflected by Eq. (15): the equation couples the observ- 
able quantities T c , Uq and £, a b- We predict for the mul- 
tilayer Hg-cuprates, Hgl212 (T c = 126 K) and Hgl223 
(T c = 135 K) the condensation energies 146 and 117 
/LteV/u.c. (£ a b « 4ao [55]), respectively using Eq. (15). 

Furthermore, we find correlation between the conden- 
sation energy Uq and ./V: the larger Uo is connected with 
smaller N (see also FIG 5). The stronger localization of 
the pair- condensate wave function seems to lead to larger 
IL condensation energy and hence to larger T c . This is 
again an unexpected result, since the stronger localiza- 
tion of the Cooper pairs should increase the Coulomb 
self-repulsion of the condensate and hence should sup- 
press T c . Note, however, that within a liquid-crystal- like 
COS the self-repulsion problem is not crucial. The hole- 
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2 4 6 

bilayer condensation energy 
FIG. 4. The critical temperature (K) at optimal doping as 
a function of the bilayer condensation energy (2(n+ l)N 2 Uo, 
meV) using Eq. (15). The real space period N is directly 
related to £ a j, via Eq. (9). The straight line is a linear fit to 
the data. The slope of the linear fit is ks which is a strong 
evidence of Eq. (15). 



anti-hole Coulomb interactions are attractive, although 
cancelled by the intra-holc and intra-anti-hole repulsions 
(intra CuOi site repulsion, ). 

Interestingly the most localized coherence area (~ 4ao) 
is provided by YBCO which is the less anisotropic mate- 
rial among HTSC cuprates (the resisitivity ratio p a b/ Pc ~ 
200 [3]). The more anistropic Hgl201 gives us some- 
what weaker localization of the condensate wave-pocket 
for similar T c . The comparison between various materi- 
als , however, is much more complicated. In general we 
can say that T c is a function of the following parameters 
considered in this study: N, p, the IL distance d and e c . 
In section VII we further analyse this complex behaviour 
of T c focusing on the calculated IL dielectric constant e c . 
The estimated xi a b values given in Tables I and II of the 
bilayer COS are in close agreement with the experimen- 
tal coherence areas t^ a b which supports the validity of our 
basic Eq. (15). 



200 



150 



■S 100 



T3 

o 



4 5 6 7 
coherence length 
FIG. 5. The measured condensation energy (given in Table 
I., /xeV/u.c.) as a function of the calculated coherence length 
(a ) using Eq. (16). 
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VII. THE CRITICAL TEMPERATURE AND THE 
C-AXIS DIELECTRIC CONSTANT 



TABLE II. The calculated dielectric constant e c using 



Eq. (18) in various c 


:uprates 


as a func 


tion of 


the coherence 


length of t 


he chai 


'ge order 


ed state. 








d(A) 


Tc(K) 


j- vn.lv ( \ 

Cb (oo) 


e c 


e c P 


CaCu02' ! 


4.64 


110 


4 


80.8 




CaCu02 


3.19 


89 


7 


83.5 




LSCO 


6.65 


39 


5 


27.9 


23 ±3, 13.5 C 








6 


27.3 










7 


11.3 




Hgl201 


9.5 


95 


5 


33.5 


34 d 


Hg- (10 GPa) 


~ 8.5 


~ 105. 


4 


27.6 




Hg- (20 GPa) 


~ 8.2 


~ 120. 


4 


25.0 




T12201 


11.6 


85 


3 


8.7 


11.3 e 








4 


26.8 










5 


13.0 










6 


8.0 




B12201 


12.2 


20 


3 


32.6 


~ 40 s 








4 


110.8 










5 


10.1 




YBCO 


8.5 h 


93 


3 


19.4 


34 / ,23.6 j 








4 


31.1 





where (,^1 c is the estimated in-plane coherence length given in 
ao ~ 3.9A. The estimated £^ !c is that value which is associ- 
ated with that calculated e c which is in good agreemant with 
the experimental e c . the bold faced values are those which 
account the best for comparison with experiment and are in 
accordance with the results of Table I. d is the CuOi plane 
to plane inter- layer distance in A [13], T c is the experimental 
critical temperature. e c is from Eq. (18). b Cao.3Sro.7Cu02 
[66], c [61], or from reflectivity measurements using Eq. (19), 
cup ~ 55cm -1 [81], A c ~ 3/im [67], d experimental e% xp values 
are taken from Am Inst, of Phys. Handbook, McGraw-Hill, 
Ed. D, E. Gray (1982), the e c value of the ionic-background 
(Hgl201: BaO), The pressure dependent T c values and IL 
distance date are taken from [71,76], e from [18], 9 from [35], 
h for YBCO the reduced IL d = 8.5A is used instead of the 
c-axis lattice constant (inter-bilayer block distance), f from 
[70], ' from reflectivity measurements: co p w 60cm -1 [81], 
A c « 0.9/im [37], e c can be deduced using Eq. (19). 



The expression Eq. (15) leads to the very simple for- 
mula for the critical temperature 



T C (N, dq,d, e c ) 



47re e c fc B 



2 N, 2N 2 

EEE 

n=l m— 2 ij 



(n) (m) 
(n,m) 



(17) 



where iVj is the number of layers along the c-axis. When 
Ni — > oo, bulk T c is calculated. T c can also be calculated 
for thin films when N, is finite and e r can also be derived 



2 Ni 2N 2 „(")_("») 

EEE 3 



47Te k B T c — ■ — — r) 

n—l m—2 tj ij 



(n,m) 



(18) 



where a c-axis average of e c is computed when Ni — > oo. 

The calculation of the c-axis dielectric constants e c 
might provide further evidences for Eq. (15) when com- 
pared with the measured values [35,61,70]. In Table II 
we have calculated the static dielectric function e c using 
Eq. (18) and compared with the experimental impedance 
measurements [64]. e c can also be extracted from the c- 
axis optical measurements using the relation [18,48] 



LOpXc 



(19) 



p auu ss c are the speed of light, c-axis plasma 



where c, iv p and A, 
frequency and the c-axis penetration depth. The general 
conclusion can be drawn on the basis of the data depicted 
in Table II that the e c values obtained from c-axis optical 
measurements are somewhat smaller then the values ob- 
tained from impedance measurements [64]. The reason 
for this is not clear. Our calculations support the opti- 
cal measurements vs. the impedance experiments if e c is 
calculated at the experimental coherence length. 

For the prototypical cuprate LSCO we get the value 
of e c = 27.3 which is comparable with the experimental 
value of 23 [61] using N = 7 which corresponds to 6ao ~ 
23.4A (£ a6 w 20-30A [79]). The overall good agreement 
of the calculated e c with the measurements is due to our 
finding that the IL charging energy is surprisingly in the 
range of ksT c when the real space period N (coherence 
length) is reasonably chosen. 

Not useless to note again the correlation in Table I. 
between £ b and Uq. There seems to be a correlation 
between the real-space localization of the Cooper wave- 
function and the SC energy gain Uq. 

The real-space "shrinking" of the Cooper wave- 
function must directly be related to the enhancement 
of pairing. Within our model the primary source of in- 
creased pairing is the IL energy gain (Coulomb induced 
pairing) . 

Optimally doped YBCO/PBCO with decoupled lay- 
ers shows no HTSC (T c w 20 K) indicating that opti- 
mal carrier concentration within isolated layers leads to 
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unit cells 

FIG. 6. The coherence length (A) as a function of the 
YBCO thickness (in unit cells) at fixed thickness of the barrier 
(insulating) layer PBCO (~ 48 and 96A) reproduced from 
Table II. of ref. [73]. 




FIG. 7. The critical temperature T c (K, Eqs. (15)-(18)) vs. 
the number of unit cells along the c-axis in YBCO using the 
4x4 model. Circles and squares correspond to the calculated 
and experimental values [5]. 



probably BCS superconductivity with low critical tem- 
perature. Setting in IL coupling (decreasing the thick- 
ness of the PBCO phase) T c is enhanced due to the in- 
crease of IL coupling [3,5]. Further experimental studies 
on artificially layered materials, such as the measurement 
of t; a b as a function of the thickness of the insulating 
phase PBCO should explain the importance of IL cou- 
pling in cuprates. The considerable increase in £ a ;, as a 
function of IL decoupling (increase in PBCO thickness 
or decrease in YBCO thickness) would be a strong evi- 
dence for the model presented in sections III.-V. Look- 
ing for such an article in the literature an interesting 
publication is found which reveals our expectation [73]. 
The extensive measurements of the upper critical field 
in (Y BCO I PBCO) n superlattices indeed results in the 
increase of ^ a b as a function of decreased thickness of 
YBCO layers. £ b varies from 4 — 5a up to ~ 8ao when 
the thickness of the YBCO phase is decreased from 10 
unit cells to 3. Unfortunatelly this article contains no 
measurements for 1 — 2 unit cell thick samples. In FIG 6 
we give the measured coherence length in YBCO/PBCO 
getting the data from ref. [73]. We predict the rapid 
increase of £ a b for dielectrically isolated single layers of 
YBCO for one or two-unit-cell thick samples which is 
the reminiscent of purely BCS features. These striking 
results represent a strong evidence for the correlation be- 
tween IL coupling and the coherence length which is 
predicted by the model presented here. Strinkingly, the 
contraction of the Cooper wave-function (the strength- 
ening of pairing) is strongly coupled to the enhancement 
of IL Coulomb coupling. 

VIII. MUTILAYER AND PRESSURE EFFECTS 

Multilayer and pressure effects on T c can also be dis- 
cussed in terms of in-plane and out-of-plane localization 
of the hole charge. In multilayer systems E\ L is com- 
posed of intra- and inter-block contributions. We use 



the notation block for the multilayer parts CaC^O^, 
Ca^Cu-^O^, etc. Eq. (14) can then be modified for mul- 
tilayer cuprates as follows, 

U « El L ^ ntra + E l c L * nter , (20) 

where 

E c L ' mtra = (I — l)E I c L {d int ra^intra), (21) 

1 denotes the number of layers. Therefore, intra-block 
charging energy further enhances T c on top of the corre- 
sponding single-layer inter-block value. YBCO is a pecu- 
liar example of cuprates in which HTSC is purely coming 
from inter-block coupling as it was explained on the ba- 
sis oiY BaiCu-sOi-?, I PrBa^CuzOi superlattice struc- 
tures [3]. In YBCO the YiCuOi)^ 1 bilayer alone does 
not show HTSC when isolated from each other in artifi- 
cially layered materials (T c w 20X). The dependence of 
T c in YBCO on the number of unit cells along the c-axis 
is calculated using Eq. (17) and the results are depicted 
on FIG 7. The experimental points [5] are relatively well 
reproduced indicating that ~ 5 unit cell thick thin film is 
already the reminiscent of the bulk properties of cuprates. 

It is worth mentioning the intercalation experiment on 
Bi2212 [4]. Intercalation of I or organic molecules into the 
(BiO)2(SrO)2 layers of Bi2212 results in no significant 
change in T c . Basically inter-layer intercalation is intro- 
duced to reduce interlayer coupling (test of inter- layer 
theory, ILT) and to enhance anisotropy in SC properties. 
Using film deposition technique, ultrathin films of Bi-Sr- 
Ca-Cu-0 (Bi2212) have been synthesized [80]. The few- 
unit-cell-thick samples show HTSC similar to that of the 
bulk material independently of the film thickness. These 
results present great challenge to ILT [10] and support 
low-dimensional SC theories. 

In the light of the N x N model, however, it is not 
surprising that in Bi2212 the nearly isolated bilayers are 
are superconductors. The single layer material Bi2201 
(Bi 2 Sr 2 CuO & ) gives very low T c (» 20K) [69], which 



10 



indicates that the dielectrics {BiO)2{SrO)2 strongly re- 
duces IL coupling, indeed the estimated e c is quite large 
(Table II) which is attributed to the weakly interact- 
ing (BiO)2 bilaycr (the BiO-BiO distance is <~ 3. 7 A 
[4]). Dielectric constant measurements [59] and the ob- 
tained small plasma frequency using c-axis optical mea- 
surements (u w bcm^ 1 ) [81] in Bi-compounds also pro- 
vide relatively large dielectric constants in accordance 
with our calculations. The extremely large conduction 
anisotropy 7 ~ 10 6 found in Bi-based cuprates also 
supports these findings [3]. Nearest-neighbouring Cu-0 
planes in Bi-compounds are nearly insulated [70] . There- 
fore, in the multilayer Bi-compounds bilayer and trilayer 
blocks are responsible for the high-T c , inter-block cou- 
pling is negligible. The multilayer Bi-compounds pro- 
vide then an example for pure intra-block HTSC. In these 
materials the multilayer-blocks are dielectrically isolated. 
However, in the most of the cuprates T c is enhanced both 
via intra- and inter-block effects. An important conse- 
quence of the bilayer N x N model is that isolated sin- 
gle layers do not show HTSC, coupling to next-nearest 
Cu0 2 plane is essential in HTSC. It is possible then to es- 
timate e c for the bilayer block, since E I c L ' mter w in the 
Bi-compounds in Eq. (20). In the hypothetical bilayer 
system in Table II {Ca(Cu0 2 )2, d = 3.191, T c = 89K), 
which is the building block of bilayer cuprates, we find 
the value of N = 8 accounting for realistic e c . N = 8 
is in accordance with the measured relatively "large" co- 
herence length of £afc w 27 A ((N - l)a = 7a w 27.31) 
[77]. 

Also, upon pressure (p) the inter-layer spacing de- 
creases, which increases E^ L without the increase of the 
hole content. Saturation of T c (p) is reached simply when 
increasing IL sterical repulsion starts to destabilize the 
system. In systems, such as YBCO or LSCO, negative 
or no pressure dependece of T c is found [3] due to the 
short Cu-0 apical distance (dcuO ~ 2 A A) which leads 
to already steric repulsion at ambient pressure and to 
the weakening of HTSC. In these systems the net gain 
in charging energy is not enough to overcome steric re- 
pulsions at high pressures. We give the calculated e c 
values in Hgl201 for high pressures in Table I. e c docs 
not depend on the pressure and therefore the pressure 
dependence of T c can be given by Eq. (17). 

With these results, we are now in a position to reach 
the conclusion that the N x N model can readily ac- 
count for at least certain physical properties of multilayer 
cuprates, such as pressure, doping and multilayer depen- 
dence of T c . Furthermore it is possible to make some es- 
timations on the upper limit of T c using Eq. (17) for T c . 
Assuming relatively small e c s=a 10 and short inter-layer 
spacing d w 7.0l we get the value of T c w 333K for a 
strongly localized electron pair with 5x5 coherence area. 
Of course, we have no clear cut knowledge at this moment 
on how the pair condensate wave-function spreads upon 
varying e c and d. Effective localization of the Cooper- 
pair wave-pocket could lead, however, to room tempera- 
ture superconductivity under proper structural and dielec- 



tric conditions if the model presented above is applicable. 

According to the stripe scenario, the charge ordered 
state of the coherence area can also be seen in many 
cuprates as one-dimensional stripe order or charge den- 
sity waves [43]. The striped antiferromagnetic order 
found in LSCO by magnetic neutron scattering experi- 
ments [42] also implies a spin-density of periodicity <~ 8a 
a reminiscent of the coherence area. The incommensurate 
"checkerboard" patterns seen with a spatial periodicity 
of <~ 8a in the vortex core of Bi2212 obtained by scan- 
ning tunneling microscopy [41] is also consistent with our 
N x N hole-anti-hole charge ordered state where N = 8 
to 9 in BSSCO. 



IX. CONCLUSION 

In this paper we studied the pair condensation and 
confinement of the hole-content on a CuC>2 superlatticc 
layer as a function of inter-layer distance and dielectric 
permittivity of the charge reservoir. 

• The assumption of a 2D 3D quantum phase 
transition of the hole-content at T c in HTSC mate- 
rials is thought to be an important general feature 
of pair-condensation and is supported by c-axis op- 
tical measurements and by first-principles calcula- 
tions. Our proposal is that the c-axis charge dy- 
namics of the hole-content contributes significantly 
to the condensation energy below T c . We find that 
the inter-layer capacitance is temperature depen- 
dent in cuprates and therefore the drop of the c- 
axis dielectric constant can be seen below T c . This 
is what leads to then the stabilization of the super- 
conducting state vs. the normal state. 

• We have found that a pair condensate can be 
distributed on a N x N square lattice layer in 
such a way that the lattice sites are filled by 
q = ±0.16e condensed charge alternatively depend- 
ing on the hole(+g)-anti-hole(— q) charge separa- 
tion dq. In this way a charge ordered state of 
the pair-condensate occurs with a "checkerboard" 
like pattern seen recently by experiment [41]. The 
phase separation of hole-electron pairs (hole-anti- 
hole pairs) in this model is stabilized electrostat- 
ically. The maximum charge separation is dq w 
0.32e, if the optimal hole content p a w 0.16e. 

• In the adjacent layers the electron-hole pairs are 
distributed in such a way that complementer charge 
pattern occurs (inter-layer charge asymmetry) in 
order to maximize the inter-layer charging energy 
A hole in one plane is neighbouring with an anti- 
hole in the other (see FIG 2.). In this way we 
derived a charge ordered N x N correlated bi- 
layer supercell model of the superconducting state 
with inter-layer charge antisymmetry which di- 
rectly leads to inter-layer Coulomb energy gain in 
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the superconducting state. The IL charge asymme- 
try could directly be tested experimentally using 
e.g. two-cell-thick thin SC films, and the "checker- 
board" STM images of both sides of the thin film 
could be measured. A complementer "checker- 
board" charge pattern would reveal our finding. 

• The superlattice nature of the pair condensate is di- 
rectly related to the smallest size of the condensate 
wave-pocket which is remarkably comparable with 
the measured in-plane coherence length of £ a & = 
10 — 20A (4ao s=s 15. 6A) of single cuprates, where 
ao ~ 3.9A is the in-plane lattice constant. The 
coherence area of the Cooper-pair wave- function in 
cuprates is strongly localized, which is due to inter- 
layer charging effects. 

• The bilayer model with 4e boson charge naturally 
implies the mass enhancement of m* w 4m e in ac- 
cordance with measurements. The calculated inter- 
layer charging energy is in the range of the ex- 
perimental condensation energy for the bilayer-hole 
system. 

• The static dielectric constant e c is calculated for a 
couple of cuprates and compared with the available 
experimental measurements. The general agree- 
ment is quite good indicating that the pure inter- 
layer electrostatic model leads to proper description 
of the static dielectric response of these layered ma- 
terials. 

• The basic microscopic mechanism of HTSC is to 
be understood within the BCS-Eliashbcrg the- 
ory. The limiting critical temperature for BCS- 
type superconductor is around 20K as it was 
found for cuprates with nearly isolated CuOi layers 
(BSCCO, YBCO/PBCO superlattices etc.). The 
detailed study of this model showed that the inter- 
layer charging energy is proportional to the ther- 
mal motion at T c , if the c-axis dielectric constant 
e c and the coherence area is appropriately chosen. 
The correlation between the coherence area £ & and 
inter-layer coupling is predicted by our model. The 
stronger IL coupling leads to smaller £ a {,. This re- 
lation is evidenced by the measurements of £ a t, in 
YBCO/PBCO films with varying YBCO thickness. 

If the physical picture derived from our model is cor- 
rect, it should be a guide for further experimental studies 
aiming to improve SC in cuprates or in other materi- 
als. This can be done by tuning the IL distance and e c 
(increasing the polarizability of the dielectric, hence de- 
creasing e c ) in these materials. The application of this 
model to other class of HTSC materials, such as fullerides 
or MgBi is expected to be also effective. 
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